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In this paper we consider a family of dynamical systems ( "arabesques" ) defined as closed chains 
of 2-element negative circuits. An n-dimensional arabesque system has n 2-element circuits, 
but in addition, it displays by construction, two n-element circuits which are both positive vs 
one positive, one negative, depending on the parity (even vs odd) of the dimensionality n. In 
view of the absence of diagonal terms in their Jacobian matrices, all these dynamical systems 
are conservative and consequently, they cannot display any attractor. First, we analyze a linear 
variant of them which we call "arabesque 0" or for short "AO" . For increasing dimensionalities, 
the trajectories are increasingly complex open tori. Next, we inserted a single, cubic, nonlinearity 
that does not affect the signs of its circuits ("arabesque 1", for short "Al"). These systems have 
three steady states, whatever the dimensionality, in agreement with the order of the nonlinearity. 
All three are unstable, as there cannot be any attractor in their state-space. The 3D variant 
("A1_3D") has been analyzed in some detail and found to display a complex mixed set of quasi- 
periodic and chaotic trajectories. Inserting n cubic nonlinearities (one per equation) in the 
same way as above, generates "A2_nD" systems. System A2_3D behaves essentially as Al_3D, 
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in agreement witli tlie fact tliat tlie signs of the circuits remain identical. System A2_4D, as 
well as other arabesque systems with even dimensionality, has two positive n-circuits and more 
(nine) steady states. It would deserve additional attention. We also investigate and compare the 
complex dynamics of this family of systems in terms of their symmetries and a connection with 
Hamiltonian dynamics is briefly presented. 

Keywords: Jacobian circuits and frontiers; Conservative systems; arabesque systems, Lyapunov 
exponents, SALI 



1. Introduction 



The fact that negative retroaction is involved in periodicity and hor neostasis was a lready clear in the works 
by Claude Bernard, who called it "elasticity", and W. B. Cannon Cannonl . 1932], who invented the term 
"homeostasis" . The notion that positive retroaction is involved in multistationarity (the occurrence of more 
than one steady state, either stable or not) is much more recent. A basis for the present understanding 
of these "nontrivial" processes is the recognition that, among the terms of the Jacobian matrix (or more 
generally in the graph of interactions) , only those terms that belong to a circuit are present, by construction, 
in the characteristic equation of the system, and hence, only these terms are involved in the kind of the 
steady states. By the term circuits we refer to those sets of terms of th e Jacobian matrix of the d ynamical 
system whose row and colu mn indices are in circular pe r mutation (see Eisenfeld &: de Lisil . Il985| ] who use 
a different terminology and Thomas fc Kaufman . 2001 ], (Thomad . 1999]). Circuits are positive or negative 
according to the sign of the product of their terms. The functions of positive and negative circuits are not 
only distinct, but contrasting. 

It has become clear that the presence of a positive circuit (in the Jacobian matrix or in the grap h 
of interactions of a system) is a necessary co nditio n for rnultista tionarity (conjecture: Thomas . 198ll |. 
formal demonstrations from (Plahte et a l. ll99.^ and [Snoussil . Il998l ] , for systems with constant sign partial 
derivatives, to [Sou 3,|2003] for the general case). It is also widely recognized that the presence of a 2- or 
more -elements negative circuit is a necessary condition for a pers istent periodici ty (conjecture by [Thomas . 
1981 ]. and demonstrated, so far limited to discrete systems by [Richardl . 201 ll ]). Concerning the number 
and nature o f attractors in r elation to circuits in systems described by differe ntial equatio i is, one can see 
for exa m ple i Tsypkin. 1962 : Richard et"al . 2011], and for Boolean networks (Coled . 19851 : Aracena et al 
2004A| Ib1: IPCTnongeo t et ail2012Af T 

In this paper, we investigate the behaviour of conservative systems (coined "arabesques") whose simple 
form can be described, through the structure of their Jacobian matrix, as closed chains of three or more 
(depending on the dimensionality n) 2-element negative circuits. An n dimensional arabesque system has n 
negative 2-element circuits, but in addition it displays by construction, two n-element circuits both positive 
vs one positive and one negative, depending on the parity (even or odd) of the dimensionality n. However, 
they can never be both negative. 

It is important to emphasize that there are no diagonal terms in their Jacobian matrices (equivalently, 
there are no 1-element circuits) and hence these systems are conservative, in the sense that their state-space 
volume is conserved (V^;^ = O), where x is a shorthand for the n dimensional vector of states {xfc}^^-^. 
In view of the absence of any diagonal element in their Jacobian matrix, one does not expect the presence 
of any attractor. 



Xk-l 

z, . . . 



These systems, of which the linear variant is briefly described in Thomas &: Rosslerl . 120061 ]. are built 
on the basis of 2-element negative circuits between variables x^+i and Xk-i as follows: = 2:^+1 - 
where k = 1, 2, 3 ... n follow a circular permutation. For convenience, we will call the variables, x, y 
rather than xi, X2, 2:3, . . . etc. 

Matrix J in Fig. [T] permits to read directly the circuits of the 3-variable linear system: it has three 
2-element negative circuits {{J12, «^2i}i {Jiz-, J'ii}^ {<^23,>^32}} and two 3-element circuits, one positive 
{{■^12, <^23, «^3i}} and one negative {{J13, J2I) "^32}}- More generally, whatever n, this construction 
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Fig. 1. The equations (a), the Jacobian matrix J (b) and the graph (c) of the linear variant for n ■ 
negative interactions are drawn in green and red respectively. 



3. In (c) tlie positive and 



implies, in addition to the presence of n 2-element negative circuits, the presence of two n-element circuits, 
either one positive and one negative, or both positive, depending whether n is odd or even. 

We would like to mention that we use the following convention for the signs of the eigenvalues of J. 
The signs of real eigenvalues can be positive, negative or zero according to the case at hand. Regarding the 
pairs of conjugated complex eigenvalues, we list behind a "/" the signs of their real part. For example, in a 
3D system, the symbol {+H — } means that all three eigenvalues are real and that two of them are positive 

(a saddle-node), while {+/ } means that there is a real positive eigenvalue and a pair of conjugated 

eigenvalues whose real part is negative (one of the two types of saddle- foci) and finally, the {-|-/0 0} means 
that there is a positive real eigenvalue and that the pair of conjugated complex eigenvalues are purely 
imaginary, i.e. their real parts are zero. 

We now pass to the next section where we present and discuss methodologies that enable one to 
characterize efficiently and fast the dynamical nature of a given orbit (i.e. discriminate between chaotic 
and regular motion) of a dynamical system, such as those we study herein. 



2. Methods for discriminating chaotic from regular motion 

In the theory of dynamical systems, the discrimination between regular and chaotic motion is of crucial 
importance since in ordered domains we have predictability while in chaotic regions, after a time period, we 
are unable to predict the evolution of the system due to the sensitive dependence of it on neighboring initial 
conditions, a conditio sine qua non signature of chaos. This is a characteristic and important feature found 
in many applied sciences, such as Physics, Chemistry, Biology and Engin eering etc., where the dynamics is 
often dissipative and the motion occurs on a strange attractor ^Sparrowl . Il982l : iLichtenberg &: Lieberman , 
1983l |. 



However, the ability to distinguish between order and chaos becomes much more difficult when the 
motion is conservative and the number of equations of motion of the system (or its dimensionality) is 
increasing, mainly since we cannot visualize their state or phase-space. It naturally follows that we need 
methods that can help us decide accurately, fast and efficiently about the chaotic or ordered character 
of a given orbit. Many methods have been developed over the years aiming to tackle this problem. The 
first of all was the well known method of the Poincare surface of section (PSS) which takes advantage of 
the intersections of an orbit with an nD "surface" (n being here the dimension of the "surface") that is 
transversal to the flow generated by the system. It can be used mainly in conservative dynamical systems 
with a small number of equations of motion, such as for example the 3D version of the arabesque systems 
we consider in this paper. For example, in 2 degrees of freedom (DOF) Hamiltonian systems where the 
phase-space is 4D, the energy integral can be used to lower the dimension to 3, thus rendering the PSS 
into a 2D surface (i.e. a plane). However, the results of this method are difficult to interpret in systems 
with more than 2 DOF or in many dimensions in general. 

There have also been developed a great variety of methods over the past years to face this subtle 
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problem, however we prefer to stick for the present study to an efficient methodology due to its su ccess to 
successfully discriminate between ordered and chaotic motion. In parti cular, in a series of papers Skokod . 
2OOII : ISkokos et aJ . l2004l : lAntonopoulos et al I2OO6I : ISkokos et all . l2003l |. the authors have used a fast, ac- 



curate and easy to compute method to distinguish regular from chaotic motion in conservative dynamical 
systems called the "Smaller Alignment Index" (SALI) method. Here, we will briefly recall its deflnition 
in subsection 12.21 and show its effectiveness in distinguishing regular from chaotic motion in the class of 
arabesque systems. Before that however, in the next subsection 12. 11 we prefer to make a short introduction 
to the theory of Lyapunov exponents that we will utilize later in the paper, since it is a standard and well 
accepted methodology for tackling the same kind of problems. 



2.1. The Lyapunov exponents 

As we have already pointed out, knowing whether the orbits of a dynamical system are ordered or chaotic 
is fundamental for understanding its behavior. In the dissipative case, this distinction is easily made as 
both types of motion are finally attracting. In conservative systems, however, the distinction between order 
and chaos is often a delicate issue, for example when the chaotic and/or ordered regions are small-sized and 
interwined, especially in multi-dimensional systems where it is not tractable to visualize their dynamics. 
For all these, it becomes crucial to make use of quantities that can decide if an orbit is ordered or chaotic, 
independently of the dimension of the state or phase-space and/or our ability to inspect visually the 
dynamics. 

The well-known and commonly used traditional method for this purpose is the evaluation of the 
maximal Lyapunov Characteristic Exponent (LCE) ci. If cji > the orbit is characte rized as chaotic . 
The theory of Lyapunov exponents was applied to characterize chaotic orbits by Oseledec [Oseledec . 1968l |. 
while the connection between Lya punov exponents and exponential diver gence of nearby orbits was given in 
Benettin etal . Il97fil : IPesinl . Il977t | . Benettin et al. |Benettin et all . I1980A|| studied theoreticall y the problem 



of the computation of all LCEs (Lyapunov spectrum) and proposed in Benettin et al 1980B| | an algorithm 
for their numerical computation. In particular, cJi is computed as the limit for t — )• 00 of the quantity: 



Li{t) = - In 



\w 



fji = lim Li{t) 



(1) 



where w{Q) and w{t) are deviation vectors from the orbit that one wants to characterize, at times t = and 
t>Q respectively. It has been shown that the above limit is finite, indepen dent of the choice of the metric 
for the state-space and converges to ai for almost all initial vectors t(J(0) [Oseledec . 19681 : Benettin et al . 
1980AI |B|]. Similarly, all other LC Es crp, cr.s etc. are comp uted as the hmits for t — )- 00 of some appropriate 
quantities, L2{t), L^{t) etc. (see Benettin et al ll980B | for more details). We note that, throughout the 
paper, whenever we need to compu te the values of the LCEs, we apply the algorithm proposed by Benettin 
et al. [Benettin et "dl . Il97fil . ll98nAl | . 

The time evolution of the deviation vectors is given by solving the variational equations (see Subsec. 
12. 2p . In general, for almost any choice of a sufficiently small initial deviation vector t(;(0), the limit for 
t —7- 00 of Eq. ([1]) gives the same a\. For reasons of avoiding numerical overflows, since the exponential 
growth of w{t) occurs for short time intervals, one has to stop its evolution after some time Ti, and record 
the computed Li(Ti), normalize vector w{t) and repeat the calculation for the next time interval T2, etc. 
obtaining flnally ai as an average over many time intervals Tj, i = 1, . . . , A'" as: 



1 ^ 



Since ai is influenced by the evolution of ty(0), the time needed for Li{T) (or Li(Tj)) to converge is not 
known a priori and may become extremely long. This makes it often difficult to infer whether o"i finally 
tends to a positive value (which signifies a chaotic behaviour) or converges to zero (which signifies a regular 
or ordered motion). 

In recent years, several methods have been introduced, which try to avoid this problem by studying 
the evolution of deviation vectors. In the present paper, as we have already mentioned, we will use the 
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method of the SALI |Skokod . [200ll | that we will describe in the foll owing subsection l2.2[ It has been applied 
successfully to several mu l ti-dimensional map s Skokos et all . 2002 ] and Hamiltonian systems Skokos et a! 
20041 : lAntonoDoulos et al l2006l : ISkokos et aJ . l2003l | and was found to converge rapidly to zero for chaotic 
orbits, while exhibits small fluctuations around non-zero values for ordered orbits. It is exactly this different 
behavior of the SALI which makes it an ideal indicator of chaoticity: Unlike the maximal LCE ai, as we 
will see in the next subsection, SALI does not start at every time step a new calculation of the deviation 
vectors, but takes into account information about their convergence on the unstable manifold from all 
previous times. 



2.2. The Smaller Alignment Index 

Let us consider the nD state-space of an autonomous conservative dynamical system: 



dx{t) 
dt 



(2) 



where x{t) = {xi{t),X2{t), . . . ,Xn{t)) and ^{x{t)) = (<I>i(x(t)), $2(^(^)), • • ■ ,^n{x{t))). The time evolution 
of an orbit associated with the initial state x(to) at initial time to is defined as the solution of the system 
of n first order differential equations (ODEs) ([2]). 

In order to define SALI we need to introduce the notion of the variational equations. These equations 
are the corresponding linearized equations of the system of ODEs ([2]) (i.e. given by the Jacobian matrix J 
of ^{x{t))), about a reference orbit x{t) and are defined by the relation: 

dwi{t) 



dt 



J{x{t))wi{t), i = l,2 



(3) 



where J{x(t)) is the Jacobian matrix of the right hand side of the system of ODEs ([2]) calculated about the 
orbit x{t) = X2(t), . . . ,x„(t)). wi{t) and W2(t) are known as deviation vectors because they describe 

small deviations of neighbor orbits from x{t) (see Eq. dZ])). We define the two Alignment Indices (ALI) as: 



and 



ALI_(t) 



ALI+(t) 



Wi{t) W2{t) 



\W2 



+ 



wim \\w2{t)\\ 



where || • || denotes the usual Euclidean norm of ^ 
SALI is then defined as: 



SALI(t) = min(ALI_(t), ALI+(t)). 



(4) 



(5) 



(6) 



It follows that one can normalize wi{t) and iv2{'t) at each integration time step without affecting their 
inner angle. It comes out that when ALI_(t) — )• 0, Eqs. dH), ^ imply that the two deviation vectors tend 
to become collinear and parallel because their inner angle tends to and so ALI_|_(t) — >■ 2, while when 
ALl4_(t) — )• the two deviation vectors tend to become collinear and anti-parallel because their angle tends 
to TT and so ALI_(t) — )• 2. 

Thus, if the orbit under consideration is chaotic then from Eq. ([6]) we get: 

lim SALI(t) = min(0, 2) = 
because both deviation vector directions tend to c oincide with the d irection of the most unstable neighbor 



Voglis et a 1. Il999l |. This implies that the two deviation 



manifold of the variational equations as t — )• cx) 
vectors tend to become colinear and either coincide or become anti-parallel to each other. 
On the other hand, if the orbit is regular (e.g. quasi-periodic) then: 



lim SALI(t) / 

t— >-oo 
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and oscillates in (0, \/2]. This is because both deviation vectors tend to become tangential to the torus on 
which the orbit is confined. The reason here is that wi{t) and W2{t) have components evolving along and 
across the torus in direction s given by two inde pendent vector fields: the Hamiltonian and a second local 
quasi-integral of the motion [Skokos et a 

This distinct behaviour of SALI for regular and chaotic orbits makes it a very clear and reliable way 
to distinguish between order and chaos in conservative dynamical systems. It is evident that the choice of 
wiito) and W2{to) is arb itrary and do e s not affect the res ults of the method. This is additionally ensured 
by results presented in ISkokosl . I2OOII : I Skokos et a I 120031 ] . Thus, along this direction, we adopt and use 
through out this paper as deviation vectors: 



tui(0) = (l,0,0,...,0), 
W2{0) = (0,1,0,...,0). 



(7) 



2.3. Circuits and Frontiers: Partition of state-space into domains according to 
the nature of the eigenvalues of the Jacobian matrix 

The state-space of a dynamical system can be partitioned into domains according to the eigenvalues of 
its Jacobian matrix J: their sign if they are real, complex or purely real, etc. Since the eigenvalues in a 
domain of the state-space determine the precise nature of any steady state present in it, th is approach can 
provide a global view of the state-space of the system [Thomas fc Kaufmanl . I2OOII . [2OO530. 



The characteristic equation, det[J — A/„] = (in which is the identity matrix in n dimensions), 
whose roots are the n eigenvalues of the Jacobian matrix J, can be written equivalently, using the two 
other standard forms: 



and 



A'^ + a„_iA"-i + a„-2A"-2 + . . . + aiA + ao = (8) 



(A - Ai)(A - A2) ... (A - A„_i)(A - A„) = 0. (9) 



Coefficients a„_j, i = 1, ... n of the characteristic polynomial (which result from the expansion of the 
determinant), are, respectively the sums of all the cyclic permutations (with the appropriate signatures, 
see below), of 1, 2, . . . , n elements of J. In other words, they are the sums of the values of the circuits and 
unions of disjoint circuits involving 1, 2, ... n variables, respectively, each with a sign (-|- vs — ) that depends 
on the parity (even vs odd) of the number of circuits involved in the union of disjoint circuits considered. 

The value of a circuit is the product of its elements. For the sake of clarity, we list the elements of a 
circuit with a simple space (implicit multiplication) but use a (multiplicative) dot to separate unions of 
disjoint circuits. For example, (an • 023032) represent the union of the 1-element circuit with value an and 
the 2-element circuit with value a23a32. 

Consequently, coefficient a„_i represents the sum of the values of the 1-element circuits, i. e. the 
diagonal terms of the Jacobian matrix, in three dimensions: —(an +022 + 033); coefficient 0^-2 represents 
the sum of the values of the 2-circuits, i.e., in 3D, — (ai2a2i+a23a32+a3iai3) and of the unions of two disjoint 
1-circuits, +(aii-a22+a22'a33 + a33'aii); coefficient ari-3 represents the sums of the values of the 3-circuits, 
-(012023031 + 013021032), of the unions of disjoints 1- and 2-circuits, +(aii • 023032 + 022 • 0031 + "33 • 012021) 
and of the unions of three disjoint 1-circuits, —{an ■ 022 • 033), etc. 

The fact that the characteristic equation is entirely built from the circuits present in the Jacobian 
matrix explains the crucial role exerted by circuits in the dynamics of systems. In fact, the arabesque 
systems treated in this work have been conceived and built in terms of circuits. It may be of interest to 
remark that the above permits to write the characteristic equation "by hand" when the number of circuits 
is rather limited and the multiplicity not too high. 



^This section represents a brief actualiz ation of the more detailed pa pers mentioned. We take this opportunity to correct a 
printing mistake found in Appendix 3 of [Thomas fc Kaufmanl . |2005A| : in the symbolic calculation program "To nDimensions 
(For Frontiers Fi, F2 and only)" one should read "eqcar = Collect[Det[- Outer[D, syst,var] + IdentityMatrix[dim]]/.values, 
] ; (a sign "-" was missing) . 
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Note that in view of Eqs. ([8]), Q, the first coefficient (the sum of the diagonal elements) is, apart 
from the sign, the sum of the eigenvalues and that the last coefficient (the determinant of J) equals (— 1)" 
the product 11 of the eigenvalues. This product 11, the coefficient oq and det(J), the determinant of the 
Jacobian matrix of the system, are related by: oq = (— l)"nAi = (— 1)"" det( J). 

The part ition of phase space according to the kind of eigenvalues is realized by utilizing the concept 
of "frontier" [Thomas fc Kaufman! . 12005 Al : lThomas fc Nardoii3 . 12nn9l |. Frontiers are defined as functions 
of the invariants of the characteristic equation of the local Jacobian matrix. The detailed algebrai c deriva- 



tion by a symbolic calcul ation program can be found in Appendix 3 of [Thomas &: Kaufmanl . I2005AI : 
Thomas &: Nardond . 2009l |. It readily translates these general operations i nto definite analytic express ions 
according to the nature of each system. Among the frontiers described in Thomas fc Nardone . 20091 ] . we 
shall use herein only frontiers Fi, i*2, which account for changes of the sign of the eigenvalues of the 
Jacobian matrix J, and -F4, which deals with the transition from real to complex eigenvalues of J. 

Fi is defined as the set of points where 

n 

ao = {-irl[X^ = {-irdet[J]=0. 
1=1 

As ao equals the product of the eigenvalues, frontier Fi is the locus of the points of state-space at which 
at least one of the eigenvalues is nil. In simple cases, one of the eigenvalues changes sign as this frontier is 
crossed, which then separates domains that differ by the sign of one eigenvalue. 

F2 is the locus of the points at which there exists a pair of purely imaginary eigenvalues of J. In simple 
cases, upon crossing this frontier, a pair of complex conjugated eigenvalues switches from two complex 
values with positive real parts (symbolized as {/ + +}) to two complex eigenvalues with negative real parts 
(symbolized as {/ }). Upon frontier F2, one has then {/O 0}. 

Finally we define frontier F4 as follows: For 2D and 3D systems, this boundary was defined analytically as 
the discriminant of the characteristic equation of J, in order to discriminate real from complex eigenvalues. 
For higher dimensionalities however, the very notion of discriminant vanishes and for this reason we define 
i<4 to be the locus of the points where there are two equal eigenvalues. Whenever a pair of complex 
eigenvalues switches to a pair of real eigenvalues, we are on such a -F4 frontier. 



2.4. Conserved quantities, symmetries and a Hamiltonian connection 

Another aspect that we need to bring forth in focus is the conservative nature of all different types of 
arabesque systems considered in this paper. Thus, let us base ourselves on the linear arabesque system and 
in particular in its 3D variant (see Fig. [TJ given by: 

dx 
dy 

— = z — X. 
dt 

dz 

—r = x — y 
dt ^ 

It is straightforward to verify and actually the same holds for the nD variant of the linear class, that the 
above system conserves a kinetic-like "energy" function C, and two Hamiltonian-type functions Hi and 
H2 which are integrals of motion: 

Hi{x,y,z) = x + y + z, (10) 
H2{x,y,z) = ]^{x^ ^y"" + z^), (11) 



C{x, y,z) = ^ {{y - zf + {z- xf + {x- yf) 



(12) 
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By construction of the class of the hnear arabesque system, the kinetic-hke energy function C{x, y, z), and 
Hamiltonian-hke functions Hi and H2 reduce to zero due to its cychc permutation symmetries. Simply, 
considering their definition and substituting by the right hand side of Fig. [H we get: 

H 

and 

dHio 



dc , dc ^ dc ^ 

dC = -r— ax + -T^dy + -r—dz 
ox oy oz 







dHi 2 = — ^iT^dx H TT^dy -\ Tr^dz 



0. 



dx dy dz dt 

Additionally, the same straightforward calculations hold for the new Hamiltonian-type function: 

1 



Ha+i{x,y,z) 



(13) 



a + 1 

where a 7^ —1, for the system, which is also a member of the arabesque larger class of systems, with all 
variables participating via simple power laws: 



dx 
dt 



-y" 



x". (14) 

^ = X- - y'^ 

Still, similar statements hold for this system, as it shares the same permutation properties as the linear 
case (see Fig. [1]). Evidently, system (jl4p conserves the kinetic-like "energy" function C and functions Hi 
and Ha+i- The cyclic permutation property is responsible for the state-space conservation condition given 
by: 

dHa+i = x'^hf - + ?/"(^" - + ^"(^° - y") = 0. 

The proof is the same as for the functions Hi and C . 

We would like to emphasize on two remarks here. First, one can easily show that the more general 
system defined by: 

dx h, J, 

y'^ - z - z +y 



dt 
dy 

dt 

dz 

IE 



X — X 



(15) 



X 



y 



/ + x 



still conserves function i^i, given by Eq. ()13p . hence we see that an "energy" dissipation or conservation 
condition, is not merely due to the presence of a nonlinearity, but depends solely on the symmetry breaking 
between the nonlinearity of the variables. 

Second, the Jacobian matrix of systems defined by E ns. (fHl) and (115]) is th e same as the discrete 
Jacobian matrix of a Boolean cascade of negative 2-circuits jAracena et all . I2004BI ] . In the case of Boolean 
tangential negative circuits of any length, the occurrenc e of limit-cycles o f long p eriod, analogous to those 
observed in the continuous case has been studied in jPemongeot et all . |2012A| |. and these limit cycles 
conserve also both the discre te kinetic energy and the global frustration f unction of the network in the 
case of a Hopfield dynamics (Pemongeot et al . 2011 : Demongeot &: Waku . 2012B ]. After this remark an 
interesting further investigation naturally suggests: That is the study of the common features observed 
in both Bo olean and cont i nuous cascades of negative loops frequently observed in genetic and epigenetic 
regulations [Cinquin et all . l2002l | . 



3. Study of the arabesque systems 
3.1. Linear variant 

In this basic version (for n = 3 and its Jacobian matrix see Fig. [T|), in view of the lack of non-linearities, the 
Jacobian and consequently its eigenvalues, are constant throughout state space. As there are no diagonal 
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elements in the arabesque systems, there cannot be any attractor. Moreover, in view of the equality of 
the coefficients used (all "1" !), the determinant of the J is nil. Consequently, the system of steady state 
equations is indeterminate, and there is thus no non-degenerate steady state, but rather steady lines (of 
equations x = y = z = . . .) 

As regarding the eigenvalues of the Jacobian matrix (see Table [T]), one or two of them are nil depending 
on the parity (odd vs even) of the dimensionality n. This can be understood by the fact that in these systems 
the characteristic equation comprises only terms with odd or even exponents depending on whether the 
dimensionality is odd or even. The other eigenvalues are pairs of conjugated purely imaginary roots. For 
sufficiently large dimensionalities (5 or > 7) two or more of these pairs display non-commensurable values, 
giving as a result non-commensurable periods: trajectories are then quasi-periodic, since they fail to close 
on themselves on the tori they generate, and thus cover it densely. The apparently aberrant simplicity 
(an ellipse, like for the cases n = 3 and n = 4) of the trajectories of the 6-dimensional system can be 
understood by the fact that the two pairs of conjugated roots have the same values (±z\/3), since their 
characteristic equation is simply lambda? {\^ + 3)^ = 0. 

A general expression of the eigenvalues of the Jacobian matrix for the nD linear arabesque system is: 

/27rA;\ 

Afe = 2isin^ j , K = 0, 1, . . . , n — 1. (16) 

This result and Fig. [3] show that linear arabesque systems of an odd dimensionality display more complex 
trajectories than the adjacent even-dimensional ones. This can be understood in terms of the presence of 
an n-element negative circuit. Like in the classical case of ellipses around a center in phase space, the size 
of the tori depends on the initial state but, in contrast with chaotic attractors, they do not display sensitive 
dependence on initial states. Fig. [1] shows projections of the trajectories (from 3 to 10 dimensions) in the 
first axes for an initial state (1, 0, 2, 1, . . .) and with Xi{{)) = 1, i > A for higher dimensionalities. 

The spectra of the Jacobians for increasing dimensionalities, clearly show this strong degeneracy ob- 
served for the non-zero eigenvalues for even dimensionalities. However, for odd dimensionalities the degen- 
eracy appears only for the zero eigenvalue as in Fig. [2j 
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Fig. 3. 3D projections (on axes x, y, z) of 3D to lOD of the linear variant (as in system of Fig. [T]). First row, from n = 3 to 
n = Q and second row, from n = 7 to n = 10. 



3.2. One cubic nonlinearity (system A1_3D) 

The next step is to introduce a nonlinearity that wih have no effect on the signs of the circuits of system 
of Fig. [H This is easily realized by rendering cubic the first term y of the first equation ^ of the system 
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Steady state 


Eigenvalues of steady state 


(-1,-1,-1) 


-0.3882 . . . , 0.1941 . . . ± 2.2612 ...i 


(0,0,0) 


0.4533 . . . , -0.2266 . . . ± 1.4677 ...i 


(1,1,1) 


-0.3882 . . . , 0.1941 . . . ± 2.2612.Z 



of Fig. [TJ We thus obtain for the A1_3D system: 



dx n 



dz 

whose Jacobian matrix is: 



— = X — y 
dt 



3i/2 -V 



J = I -1 1 
1 -1 

(which we now call "arabesque 1_3D" and for short Al_3D). Consistent with the existence of positive 
circuit(s), this nonlinearity results in the occurrence of three unstable steady states, as expected from the 
absence of any attractor. 

The determinant of the Jacobian matrix is no more nil and so the system of steady state equations is no 
more indeterminate. Moreover, consistent with the simultaneous presence of a nonlinearity and of positive 
circuit (s), the system now displays multistationarity meaning that there are three colinear steady states, 
all unstable as expected by the absence of any attractor. As we have already pointed out, the state-space 
of system (I17|] can be partitioned into domains according to the nature of the eigenvalues of its Jacobian 
matrix (see Sec. 12. 3p . We recall that frontier Fi is the locus of points at which one or more eigenvalues of 
J are nil, F2 is the locus at which there is a pair of purely imaginary complex conjugated eigenvalues of J 
and frontier F4 is the locus of points where two or more eigenvalues of J are equal. For the system (I17p . 
the analytic descriptions of these frontiers take the form: 

Fi : -1 + 3y2 = 



F2 
Fa 



-1 + 3y2 = 
27(1 -3?/2)2 + 4(2 + 3y2)3 = o 



We observe that Fi and F2 coincide. These frontiers define two planes given by y = ±-^. Between them, 

the signs of the eigenvalues are of the type {+/ }. Outside them, the signs of the eigenvalues are of 

the type {— / + +}• There is thus at the level of these frontiers, a shift from + to — of a real eigenvalue 

(i.e. the effect of frontier Fi ) and a shift from to ++ of the real part of a pair of complex conjugated 

eigenvalues (i.e. the effect of frontier ^2)- There is no F^ frontier since function F4 is strictly positive for 
all real y, in agreement with the fact that the eigenvalues remain complex everywhere in the state-space 
of system p!7|) . 

In order to get a picture of the dynamics of the nonlinear arabesque (jl7p . we have performed a series 
of numerical integrations of a large ensemble of initial states and found that its trajectories are confined in 
two compact domains which are located symmetrically to the origin which is the center of symmetry. Due 
to its special symmetry properties, to each trajectory starting from an initial state (xq, yo, zq) corresponds 
an "opposite" trajectory starting from the opposite initial state {—xq, —yo, —zq), which is symmetrical with 
respect to the center of symmetry (i.e. the origin). These point-symmetric domains lie between the two 
"outer" steady states (1, 1, 1) and (—1, —1, —1). 

We show typical time evolutions of three such characteristic orbits of (jl7p in Fig. HI They are seen to 
originate, respectively, from initial states located (a) inside (i.e. {x,y,z) = (0.7,0,0)), (b) at the periphery 
(i.e. {x,y,z) = (—0.2,0.6,0)) and (c) outside (i.e. {x,y,z) = (0.3,0,0)) of the compact domains mentioned 
above. In contrast with trajectories (a) and (b), the third trajectory escapes to infinity at around t = 500, 
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Fig. 4. Time evolution of variables x, y and z of system (|17[) superimposed with different colors, for the initial state: (a) 
(0.7, 0, 0) which is located inside the compact domain, (b) (—0.2, 0.6, 0) which lies at the periphery of the compact domain 
and (c) (0.3, 0, 0) which is located outside the compact domain, between t — 29900 and t = 30000 for panels (a) and (b), and 
between t — 300 and t = 400 for panel (c). Evolution of the maximum Lyapunov exponent ai (black) and second maximum 
o"2 (gray) for: (d) the initial state of panel (a), (e) the initial state of panel (b) and (f) the initial state of panel (c). 




Fig. 5. xy plots of trajectories of initial states (0.7, 0, 0), (—0.2, 0.6, 0) and (0.3, 0, 0) (red) and of their "opposites" (—0.7, 0, 0), 
(0.2,-0.6,0) and (-0.3,0,0) (blue). 



and that is why it is shown only from t = 300 to t = 400. In panels (d) to (f) of Fig. HI we present the 
time evolution of the corresponding Lyapunov exponents ai and a2 of cases (a) to (c) respectively. We 
deduce that the orbit which is initially located inside the compact domain (see panel (a) and (d) of Fig. 
HI is ordered at least up to t = 3 • 10^ since its maximum Lyapunov exponent ai decreases towards zero 
following a power law oc j, the orbit which is initially located at the periphery of the compact domain (see 
panel (b) and (e) of Fig. H]) shows a chaotic behavior since its ai converges to a positive value at least up 
to t = 3 • 10^ and finally the orbit which is initially located outside the compact domain (see panel (c) and 
(f) of Fig. |4]) presents a transient chaotic behavior before it escapes to infinity at about t = 400. 

The xy plots of the trajectories that correspond to the plots of panels (a), (b) and (c) (red) of Fig. |4] 
and their "opposite" (blue) are shown in panels (a), (b) and (c) of Fig. [5l 

Fig. [6] shows the xy plot from t = 29900 up to t = 30000 of the trajectory starting from the initial 
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Fig. 6. xy plot of the trajectories of system (|17|l starting from initial states (—0.2,0.6,0) (red) and (—0.2,0.6,0.0000001) 
(blue), from t = 29900 to t = 30000. The time evolution of the Lyapunov exponents cri and (J2 of the initial state (—0.2, 0.6, 0) 
is given in panel (e) of Fig. |4l 



state (—0.2,0.6,0) which is located at the periphery of the compact domain, superimposed to the xy plot 
of the trajectory starting from the neighboring initial state (—0.2,0.6,0.0000001), which differs only at the 
seventh decimal digit (from t = 29900 up to t = 30000) of the z coordinate. The obvious difference in their 
time evolutions indicates an extreme sensitivity to the initial conditions, thus strongly suggesting a chaotic 
motion as we have already pointed out when discussing the results of panel (e) of Fig. [4] regarding the time 
evolution of the maximal Lyapunov exponent ai. 

An xy plot of a typical pair of "opposite" trajectories with initial states (—0.2, 0.6, 0) and (0.2, —0.6, 0) 
is shown in Fig. [6] and in Fig. [7] an xyz plot of the same trajectories. 

In order to have a thorough picture of the dynamics of the conservative (jl7p arabesque system, we 
iterated in Fig. [8] 10^ initial states located on a rectangular equally spaced grid for (x,y) G [—1.5,1.5] x 
[—1.5, 1.5] with z = up to final integration time tf = 10^. We computed their SALI(t) values following 
the methodology described in Subsec. 12.21 and focused on those that do not escape to infinity earlier than 
tf. We categorized the non-escaped orbits as follows: if SALI(tj) > 10~^ then we characterize the initial 
state as regular (quasi-periodic or periodic) and color it gray in the plot of Fig. ([8]), otherwise (i.e. if 
SALI(tj) < 10~^) as chaotic and color it black. White color corresponds to initial states that escape to 
infinity before tf. In order to apply this kind of binary filtrat ion, we used as a chaos threshold the value 
10^^ as it has been used successfully in previous publications [Skokosl . I2OOII : ISkokos et al . l2004l ] as a good 
one for the efficient disc riminat i on between regular i ty arid chaos in Hamiltonian system s, symplectic maps 
and galactic potentials (Skokod . 2001 : Skokos et al . 2004 : CapuzzoDolcetta et al . 2007 ]. 

In Fig. [SI we can see that there are two point-symmetric "kidney" -shaped big regions of regular motion 
each one bounded by its own point-symmetric thin chaotic layers. These chaotic layers cover the region 
of regular initial states and create complex structures of chains of stability islands of small size and of 
high period. Outside the two main big regular regions there are smaller stability islands of higher period 
which are surrounded by their own even thinner chaotic layers. It is interesting to note that orbits starting 
initially from the interior of the regular region (i.e. inside the "kidney" -shaped region of Fig. [8]) are not 
able to penetrate the surrounding thin chaotic layer and subsequently escape to infinity later on, and so 
this property permit us to refer to the "kidney" -shaped regions as compact regions. 

In Fig. [9] we present the corresponding PSS of points (x, y) of system (llTp for 2; = 0. In particular, 
we see the special state-space structures (i.e. foliated tori) of the two main regular gray point-symmetric 
regions of Fig. [8] which are not however visible in this figure due to the fact that we have only used a binary 
rough filtration which is able to distinguish between order and chaos only. Interestingly enough, inside each 
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Fig. 7. Plot of the trajectories of the initial states (—0.2,0.6,0) (red) and (0.2,-0.6,0) (green). The three colinear steady 
states are seen as black points. 

"kidney" -shaped region, there are two main ehiptic points (the centres of the big ehiptic regions) separated 
by a hyperbohc one. Around these elhptic points we find again chains of stabihty islands of high period and 
of relatively big size depicted as smooth elliptic closed curves. This special structure around the elliptic 
points is separated by the separatrix coming out of the unstable hyperbolic point. As we move away from 
the two main elliptic points towards the origin (point of symmetry), we can see similar in nature but even 
more complex thiner structures of chains of stability islands of smaller size and of higher period. This is so 
until we reach the complex "border" of the "kidney" -shaped region, the so-called "edge of chaos", outside 
of which there is a narrow chaotic region filled by randomly scattered points and initial states that escape 
to infinity. 

To pursue this finding further, we have first approximated the layer of all initial states (x, y, z) of 
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Fig. 8. Cross section of points (x, y) with z = of the state-space of the nonlinear arabesque system (|17p using the SALT 
method. Each point represents an initial state integrated up to = 10^. Gray points represent regular (quasi-periodic or 
periodic) initial states with SALI(fy) > 10~*, black chaotic initial states with SALI(ty) < lO"'^ while white corresponds to 
initial states that escape to infinity before tf. 

J 




Fig. 9. PSS of points {x,y) satisfying 2 = of the state-space of system (|17|) as a juxtaposition for 2 > and 2; < due to 
its point-symmetric property. 

system (jl7p which lead to chaotic behavior. This was achieved by considering initial states on an equally 
spaced cubic grid {x,y,z) £ [—1.5,1.5] x [—1.5,1.5] x [—1.5,1.5], keeping only those initial states whose 
trajectories are characterized by a value of SALI(tj) < 10^® (see Fig. [10] (a)). Then, we computed the 
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Fig. 10. (a) Volume of chaotic initial states of the nonlinear arabesque system p7|l . We have located them keeping only those 
that give SALI(t^) < 10~*. (b) Volume of regular initial states of the same arabesque system. We found them keeping only 
those that give SALI(tj) > 10~*. In this plot we considered as final integration time tf — 10^. 



fractal dimension Sarraille et al 1994 ] of this set of points that approximates the thin chaotic layer and 
found that it is approximately 2.4 ± 0.025. Such a value signifies a fractal object clearly located between 2 
and 3 dimensions. By doing the same for one of the two symmetric main volumes of regular initial states 
(see Fig. [TT|) we have found that its fractal dimension is numerically very close to a three dimensional 
object since FD3 program reports 2.98 ± 0.3 (see Fig. [10] (b)). 

The fact that most of the trajectories starting from the interior of the point-symmetric "kidney" -shaped 
regions are quasi-periodic and most of those on the edge of the "kidney" -shaped regions chaotic, implies 
that these trajectories are confined in the domains they come from and consequently do not "merge" with 
each other. One can thus expect that trajectories starting from inside the "kidney" -shaped regions should 
fill these two point-symmetric volumes and consequently these volumes should be compact. This could be 
checked by computing the corresponding PSS, that is a selection of points of trajectories (one color for 
each individual trajectory) that lie for example on the plane (x, y) with z = (see Fig. [11] (a) and (b)). 

The two panels of Fig. [11] give an idea of the complex structure of the edge and of the core of the 
"kidney" -shaped regions. The left panel of Fig. [11] shows the points of a chaotic trajectory, located in 
the edge while the right panel of Fig. [Til those of a quasi-periodic trajectory, located in the core of the 
"kidney" -shaped regions. Systems of type Al_nD with n > 3 have been studied in less detail. For odd 
dimensionalities, they behave essentially, mutatis mutandis, like Al_3D, including the number of steady 
states, which remain three, consistent with the degree of the unique nonlinearity. Systems Al_nD with 
even dimensionalities have two nil eigenvalues, like the linear systems. The determinant of the Jacobian 
matrix is nil, and consequently the steady state equations of the system are indeterminate. There are in 
fact three steady lines in 4D (i.e a line x = z ui each of the planes y,u = —1, 1, 1) and in 6D there are 
three lines x = z = u etc. This multi modality is understandable in view of the simultaneous presence of a 
nonlinearity and of positive circuits. 

It is of interest to remark that in these even-dimensional systems, it was possible to eliminate the 
positive circuits simply be permuting the signs of the elements of any of the 2-element circuits. As the 
determinant of J is no more nil, the system is no more indeterminate, and there is a single steady state, 
in view of the absence of a positive circuits. 
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Fig. 11. (a) PSS on the plane {x, y) with 2: = of five pairs of "opposite" trajectories (from t = up to t = 1000) of system 
H17p. starting from inside the "kidneys", (b) PSS on the plane {x,y) with 2: = of five pairs of "opposite" trajectories starting 
from the edge of the "kidney" -shaped regions. Each color corresponds to points of the same trajectory. 



3.3. Systems with n cubic nonlinearities (A2-3D & 
3.3.1. System A2_3D 

In this subsection, instead of inserting an extra nonlinearity to the second variable in the first equation, as 
we have done in the previous section, we apply a cubic nonlinearity to each positive term of the arabesque 
system (jl7p . This slightly alters the simplicity of the system, but, on the other hand, increases the number 
of its symmetries. In 3D, the equations read: 

dx o 

1 = .'-.. (18) 

dZ n 

— = X — y 
dt ^ 



and its Jacobian matrix is: 



32/2 -r 
J = I -1 3^2 

.3x2 _i Q 



The presence of three nonlinearities in (llSp does not increase substantially the complexity of its behavior 
and similarly to system p!7|) . there are still three steady states (— 1,— 1,— 1), (0,0,0), (1,1,1) with its 
trajectories being confined again inside two "kidney" -shaped volumes. The analytical expression of the 
new frontiers -Fi, F2 and F/^ are: 

Fi : -l + 27x2y2^2 ^0 

F2 : -l + 27x2y2^2 ^0 

F4 : 4(3x2 + 3y2 + 3^2)3 + 27(1 - 27x2^2^2)2 = q 

Frontiers -Fi = and F2 = are no more simple planes, but eight hyperbolic surfaces (see Fig. [T2]) . 

However, like system (jl7p . there are only two types of eigenvalues of J: either {+/ } or {— / + +}. This 

is why frontiers Fi = and F2 = coincide and there is no frontier F4 = as it is strictly positive for all 
X, y and z. These similarities between systems ()17p and (jlSp are not surprising, as far as one realizes that 
the signs of all circuits involved are of the same type in the two systems (see Fig. [T2]) . 
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Fig. 12. Plot of frontiers Fi and F2 of system (|18p which coincide. In the region confined by the two frontiers, the real 
eigenvalue is negative and the real parts of the complex eigenvalues are positive {— / + +} like it is for the steady state (0, 0, 0). 
Outside the eight hyperbolic surfaces, the real eigenvalue is positive and the real parts of the complex ones are negative 
{+/ }, like at the steady states (1, 1, 1) and (—1, —1, —1). 



Steady State 


Eigenvalues of steady state 


(-1,-1,-1) 


2,-1 ± 2^/2,1 


(0,0,0) 


-1,0.5± 


(1,1,1) 


2,-1 ± 2V3i 



Table [3] summarizes the eigenvalues of system (jlSp at its steady states and Fig. [12] shows frontiers Fi 
and F2 which coincide. 



3.3.2. System A24D 

The 4D case of the nonlinear arabesque system which we call for short A2_4D is given by: 



and its Jacobian matrix by: 



y 



dx 
dt 

m = 

iz — y3 



dt 

( 



u 

X 

y 

■ z 



(19) 



J 



-1 




3y2 
3z2 

-1 
-1 



-1\ 


\3x^ -1 y 

Here and in fact, in all A2 systems with an even dimensionality, both n-element circuits are positive, and, 
in agreement with this, the number of steady states increases from three in A2_3D to nine in A2_4D. 
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Steady state 


Eigenvalues of steady state 


(-1,-1,-1,-1) 


±4z, ±2 


(-1,0,-1,0) 


±2l, ±y/2l 


(-1,1,-1,1) 


±Ai, ±2 


(-1,1,-1,1) 


±2i,±^/2i 


(0,0,0,0) 


±z,±l 


(0,1,0,1) 


±2l,±^/2l 


(1,-1,1,-1) 


±4^, ±2 


(1,0,1,0) 


±2i,±^j2i 


(1,1,1,1) 


±4«, ±2 



Table [4] shows these steady states and their corresponding eigenvalues. As one can see, the steady states 
can be classified into three groups regarding the signs (and in fact the values) of the eigenvalues of J. 
Steady state (0, 0, 0) has two real, namely ±1, and two purely imaginary, namely itz, eigenvalues. The four 
"external" steady states (i.e. those whose coordinates are all ±1) also have two real, namely ±2 and two 
purely imaginary, namely ib4i, eigenvalues. The "intermediate" steady states (i.e. those whose coordinates 
are ±1 or 0) display two pairs of conjugate purely imaginary eigenvalues, namely ib2z and ^\f2i. 
The analytical expression of the frontiers of system (jl9p are: 

Fi : -1 + 9ti2y2 + 92;2^2 _ g;^^2y2^2^2^2 ^ Q 

-F2 : 

F4 : -l + 9n2y2 + 9a;222_8]^2,2y2^2^2 ^0 

As one can see, F2 is nil everywhere (i.e. there is no F2 frontier), and frontiers F\ and -F4 coincide. This 
describes the fact that there are complex eigenvalues everywhere, and that two adjacent domains of the 
state-space of system ()19p differ. At the level of the double F1-F4 frontier, one pair of real eigenvalues 
vanishes and beyond this limit it is replaced by a pair of conjugated, purely imaginary eigenvalues. 

The steady states can be visualized, for example, by sections following the three planes (y = — 1, u = 
— 1), (y = 0, u = 0) and (y = 1, u = 1). In either of these planes, there are three steady states whose x 
and z coordinates are (—1, —1), (0,0) and (1, 1). 

Figs. [T3I and [Til show respectively in plane (y = —1, u = —1) and in plane (y = 0, u = 0), (a) the 
frontiers and steady states and (b) trajectories from initial states close to steady states (these trajectories 
are also seen in (a) at the level of the relevant steady states). The situation in plane (y = 1, n = 1) has not 
been represented since it is simply symmetrical to that of plane (y = — 1, u = — 1). One sees that compact 
trajectories (in red) are generated from initial states close to the four steady states whose coordinates 
include both and ±1 (thus, the four steady states whose eigenvalues are two pairs of conjugated purely 
imaginary eigenvalues). 

4. Conclusions and discussion on open problems 

In this paper we have studied a family of dynamical systems, named "arabesques" , whose Jacobian matrices 
possess closed chains of 2-element negative circuits. In view of the absence of diagonal terms in their 
Jacobian matrices, all these dynamical systems are conservative and consequently, they do not posses 
attractors. First, we have analyzed the linear variant of this family and then, we have inserted a single, 
cubic nonlinearity in one of its equations so that the signs of its circuits are not affected. We have analyzed 
the new nonlinear system in some detail and found that it displays a complex mixed set of quasi-periodic 
and chaotic trajectories. Next, we have inserted n cubic nonlinearities, one per equation, and generated 
in particular the A2_3D arabesque system. We have shown that it has nine steady states in agreement 
with the order of the system and that it behaves essentially as Al_3D does, since they share the same 
circuit signs. We then passed to the system A2_4D, that has two positive 4-circuits and discussed about its 
dynamics in terms of its Jacobian circuits and frontiers. Finally, we have also investigated and compared 
the complex dynamics of this family of dynamical systems in terms of their symmetries and have presented 
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Fig. 13. (a): Section showing the plane (y = —1, u — —1) with its three steady states 1, — 1, — 1, — 1), (0,-1,0,1), 
(1, —1, 1, —1) (black points). The blue curves are the section of frontiers Fi and ^4. (b) From initial states in the vicinity of 
state (0, —1,0, —1), one gets closed trajectories, one of them shown in detail in the right part of the figure. 




Fig. 14. (a) Section showing the plane {y — 0, it = 0) with its three steady states (—1,0, —1,0), (0,0,0,0), (1,0, 1,0) (black 
points). The blue curves are the section of frontiers Fi and F4. (b) From initial states in the vicinity of states (—1,0, —1,0) 
and (1, 0, 1, 0) one gets regular trajectories, one of them shown in the right part of the figure. 



Fig. 15. Along a cross-section of planes {y = 0, u = 0) from {x,y) = (—1, —1) to {x,y) = (1, 1), Fi, F4 frontiers are crossed 
at two places. One pair of purely imaginary eigenvalues is present throughout the section. The other one vanishes at the level 
of the frontier, beyond which it is replaced by a pair of real eigenvalues (turquoise). 

a connection with Hamiltonian dynamics. We conclude here by mentioning some open problems which are 
related to our work and we think they deserve to be addressed. 
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The first one is the role of functions C and H^+i, i.e., kinetic energy-hke and Hamiltonian-hke functions 
in Eq. (|12p . where we see that an "energy" dissipation or conservation condition is not merely due to the 
presence of a nonlinearity, but depends solely to the symmetry breaking between the nonlinearity of the 
variables. In view of our previous expertise in Hamiltonian dynamics via SALI, Lyapunov Exponents and 
fractal dimensions, we propose as future line of research to investigate the role of circuits in classical and 
quantum mechanical Hamiltonians (like Bose-Einstein condensate versus Fermi-Pasta-Ulam Hamiltonian 
systei ns) as well as to study th e influence of a small additional potential term (e.g., due to an external 



field) [Demongeote et al 1 . 12007 1 . 



The second is related to the case of Boolean tangential negative circuits of any length where we 
have observed th e occurr ence of limit-cycles of long period, analogous to those of the continuous case 
Demongeot et al 2012A |. These limit cycles conserves both the dis crete kinetic energy and the global 



frustration function of the network, in the case of a Hopfield dynamics [Demongeot fc Waku . 2012B]. After 
this observation, an interesting further investigation suggests by itself: The study of common features 
observed in both Boolean and continuous cascades of negative loops frequently observed in genetic and 
epigenetic regulations and try to extend the present results to circuits in tersecting in two or mo re points, 
showing a smaller number of invariant orbits than in the tangential case [Demongeot et al . 2011]. 

The last one is that symmetries found in the Jacobian matrix affect the eigenvalues and frontiers and 
a comparison with symmetries observed in known chaotic Chuas circuits and symplectic maps could be 
further envisioned in order to check if their influence on the eigenspectrum is of the same nature. 
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